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An adiabatic-connection fluctuation-dissipation theorem approach based on a range separation of 
electron-electron interactions is proposed. It involves a rigorous combination of short-range density 
functional and long-range random phase approximations. This method corrects several shortcomings 
of the standard random phase approximation and it is particularly well suited for describing weakly- 
bound van der Waals systems, as demonstrated on the challenging cases of the dimers Be2 and Ne2. 
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Density functional theory (DFT) is a powerful ap- 
proach for electronic-structure calculations of molecular 
and condensed-matter systems [1]. However, one diffi- 
culty in its Kohn-Sham (KS) formulation using local den- 
sity or generalized-gradient approximations (LDA and 
GGA) is the description of non-local correlation effects, 
such as those involved in weak van der Waals complexes, 
bound by London dispersion forces [2], The adiabatic- 
connection fluctuation-dissipation theorem (ACFDT) ap- 
proach is one of the most promising ways of constructing 
highly non-local correlation functionals. This approach, 
introduced in wave function theory [3] and in DFT [4, 5] , 
consists in extracting non-local ground-state correlations 
from the linear charge density response function. 

Recently, the ACFDT approach has received renewed 
interest for implementing the random phase approxima- 
tion (RPA) or other related approximations for atoms, 
molecules and solids [6-11]. The RPA correlation energy 
is consistent with the use of the exact, self- interaction- 
free exchange energy. In spite of a number of encour- 
aging results, such as the correct description of dis- 
persion forces at large separation [12], the proper re- 
production of cohesive energies and lattice constants of 
solids [10, 13, 14] and an improved description of bond 
dissociation [6, 7, 15], several aspects of the RPA are still 
unsatisfactory. 

First, the RPA is a poor approximation to short-range 
correlations, leading to correlation energies that are far 
too negative [16]. Second, in a Gaussian localized basis, 
RPA calculations have a slow convergence with respect 
to the basis size [6] . Third, the presence of an unphysical 
maximum (bump) at medium distances in dissociation 
curves of simple diatomic molecules [6, 15] indicates an 
inherent problem which has not yet a fully clarified origin. 
Fourth, although in principle the orbitals should be cal- 
culated self-consistcntly [17], most RPA implementations 
consist of a post-KS single-iteration calculation, making 
the choice of the input orbitals sometimes critical. Last 



but not least, although the main advantage of the RPA is 
supposed to be the description of dispersion forces, rare 
gas dimer potential curves calculated from LDA or GGA 
orbitals are often qualitatively wrong, as shown later. 

The poor short-range behavior can be corrected by 
adding a GGA functional constructed from the differ- 
ence of the exact and RPA correlation energies of the 
uniform electron gas [16], but this so-called RPA+ tech- 
nique does not lead to consistent improvement [6]. One 
can go beyond RPA by including exchange-correlation 
(xc) kernels [8, 9, 18], but so far it remains imperfect, 
e.g. local xc kernels produce pair densities that diverge 
at small interparticle distances [9, 18]. 

In a similar spirit as in the work of Kohn et al. [19], we 
propose an ACFDT approach based on a range separa- 
tion of electron-electron interactions. It involves a rigor- 
ous combination of a short-range density functional with 
one of the possible long-range generalizations of the RPA. 
The method offers a solution for several of the aforemen- 
tioned difficulties and is particularly well suited for the 
description of weakly-bound van der Waals systems. 

Theory. In the range-separated multideterminant ex- 
tension of the KS scheme, an alternative approach to 
DFT (see, e.g., Ref. 20), the exact ground-state energy 
of an electronic system is expressed as 



E = min 
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dn{(*|f + V ne + + (1) 

where T is the kinetic energy operator, V ne 
is the nuclei-electron interaction operator, 
W* = (1/2) //dridr 2 ^ e (n 2 )n 2 (n,r 2 ) is a long- 
range electron-electron interaction written with 
w ee( r ) = erf (/ir)/r and the pair-density operator 
n 2 (ri,r 2 ) = ra(ri)n(r 2 ) - n(ri)£(ri - r 2 ), and E^ xc [n] 
is the corresponding ^-dependent short-range Hartree- 
cxchange-correlation (Hxc) density functional that 
Eq. (1) defines. The minimizing multideterminant wave 
function, denoted by , 5 lr , corresponds to a long-range 
interacting effective Hamiltonian and yields the exact 
density. The parameter /i in the error function controls 
the range of the separation. For fi = 0, the standard 
KS scheme is recovered: w' r e (r) vanishes, "P 1 ' reduces 
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to the non- interacting KS wave function and E^ xc [n] 
becomes the usual Hxc functional. For fi— >oo, the usual 
wave function formulation of the electronic problem is 
retrieved: w'' e (r) becomes the full Coulomb interaction, 
Eft xc [n] vanishes and 4 ,lr becomes the exact ground- 
state wave function. For intermediate values of /x, the 
interaction effects are divided between the long-range 
interacting wave function \I/ lr and the short-range den- 
sity functional £?Hxc[ n ]' an< ^ one ex P ec ts to find better 
approximations for each piece. Short-range LDA [21, 22] 
and several beyond-LDA approximations [20, 23-25] 
have been proposed for E^ xc . Here, we use an ACFDT 
approach for the long-range part of the calculation. 

In a first step, the minimization in Eq. (1) is restricted 
to single-determinant wave functions leading to a 
range-separated hybrid (RSH) scheme [26] 

E RSH = mm {(*|f + V ne + W*\$) + (2) 

which, in contrast to some range-separated KS 
schemes [27, 28], does not include long-range correla- 
tion. The minimizing determinant <&o is given by the 
self-consistent Euler-Lagrange equation 

(f + V ne + Vf^HF + VbIc) l*o> = £b|*o), (3) 

where V^ x HF is a Hartree-Fock (HF) type long- 
range Hartree-exchange (Hx) potential, Vff xc = 
J dr5Eft xc [n<s, ]/5n(r)n(r) is the short-range local Hxc 
potential and £o is the Lagrange multiplier for the nor- 
malization constraint. As usual, V^ x HF is the sum of a 

local Hartree part V H '^ and a non-local exchange part 

v x,HF- 

The RSH scheme does not yield the exact energy and 
density, even with the exact short-range functional E^ xc . 
Nevertheless, the RSH approximation can be used as a 
reference to express the exact energy as 

E = E RS k + E%, (4) 

defining the long-range correlation energy E l * , for which 
we will now give an adiabatic connection formula. For 
that, we introduce the following energy expression with 
a formal coupling constant A 

E x = min{(*|T + V ne + V$ x HF + AVF lr |*) 

where ^ is a multideterminant wave function, W lT is the 
long-range fluctuation potential operator 

W h = WH - V& iHP) (6) 

and E^ xc is the previously-defined A-indepcndcnt short- 
range Hxc functional. The minimizing wave function is 
denoted by ^\ r . For A = 1, the physical energy E — E\=i 



and density are recovered, as Eq. (5) reduces to Eq. (1). 
For A = 0, the minimizing wave function is the RSH 
determinant ^\ r =0 = $o- Note that, because the density 
at A = is not exact, the density is supposed to vary 
along this adiabatic connection. Taking the derivative 
of E\ with respect to A, noting that Ex is stationary 
with respect to and reintegrating between A = and 
A = 1 gives 

E = Ex=o + [ d\ (^\W lT \^), (7) 
Jo 

with E x =o = (*o|r+K e +^,HFl*o)+-B H r x C [«*o] =£RSH- 
(^ol^ 11 l^o)- Thus, the long-range correlation energy is 

E% = j\\{(^\W h \*\ T ) - («fo|^ h '|$o)}, (8) 
or, using a compact notation, 

E c=lJ o d\Tr[w h *P* x ], (9) 

where w h and P^ x are four-index representations of the 
fluctuation potential and the correlation contribution of 
the two-particle density matrix in a one-electron basis, 
* stands for contraction of two indices and Tr is the 
trace over the remaining two indices. The fluctuation- 
dissipation theorem is then used to express P^ T X with the 
imaginary-frequency four-point polarizability x'^(iu) cor- 
responding to the wave function (see, e.g., Ref. 3) 

1 r°° 

Pc,x = -5- / due- M+ [ Xx T (iu) - X o(i«)]+A?, (10) 

where Xo(iu) is the four-point polarizability for the RSH 
effective Hamiltonian of Eq. (3) and is the con- 
tribution coming from the variation of the one-particle 
density-matrix along the adiabatic connection. The ex- 
pression of Alj' is straightforward but it is sufficient to 
write it as A l x T = F[G^] — F[Gq] where F is a known 
functional, G 1 ^ is the two-point one-particle Green func- 
tion corresponding to the wave function and Go is 
the two-point RSH Green function. Along the adiabatic 
connection of Eq. (5) , the Green function G 1 ^ satisfies a 
self-consistent Dyson equation 

(G l x) = G 1 - A (EhJG 1 ^] - EhJCo]) ~ ^Jf A 

(11) 

where AE^ and Sj. 1 ^ are the Hartree-exchange and cor- 
relation self-energies associated with the long-range inter- 
action uigg. The long-range polarizability is given by the 
solution of the Bethe-Salpeter-type equation (see Ref. 29) 

(x 1 a)" 1 = (x1p,a)" 1 -^-/^, (12) 

where Xjp x i s an independent-particle (IP) polarizabil- 
ity whose expression is a frequency convolution of two 
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FIG. 1: (Color online) RPA, RPAx and RSH+RPAx total en- 
ergies of Ne2 at equilibrium distance (5.84 Bohr) with respect 
to the basis size. 



Green functions G A , abbreviated as x}p A = — iG A G A , 

and Xf^^XST^/SG^ and fc.x = ^c.a/^a arc lon g" 
rangc HF-type and correlation kernels. 

So far, the theory is in principle exact. In the following 
we introduce the approximation 
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which corresponds to neglecting long-range correlation in 
all one-electron properties. Indeed, from Eq. (11), this 
approximation implies that the Green function remains 
unchanged along the adiabatic connection, i.e. G A = Go 
and thus A A = 0. It also follows that /' r A = and 
Xip a = Xo, and Eq. (12) then has the structure of the 
RPA with HF exchange kernel (sometimes refers to as lin- 
ear response time-dependent Hartrce-Fock theory or full 
RPA) that we will designate here by RPAx, by opposi- 
tion to standard RPA, without exchange kernel (some- 
times also called direct RPA). In the basis of RSH spa- 
tial orbitals, for spin-restricted closed-shell systems, the 
long-range correlation energy then becomes 



E. 



iajb 



\ab) 



iajb ' 



(14) 



where ia and jb refer to excitations from occupied (i 
or j) to virtual (a or b) orbitals, (ij\w l * e \ab) are the 
two-electron integrals with long-range interaction, and 
(Pl T x )iajb are the matrix elements of the spin-singlct- 
adapted P^. The one-electron term V# x HF in Eq. (6) 
does not contribute to E 1 * because of the occupied- virtual 
structure of P^x- O m y singlet excitations contribute to 
Eq. (14), since the two-electron integrals involved van- 
ish for triplet excitations. Note that alternative (but 
inequivalent) RPAx correlation energy expressions, such 
as the plasmon formula of Ref. 3 and the closely related 
ring CCD approximation of Ref. 30, require contribu- 
tions from both singlet and triplet excitations, which 
may be problematic in systems displaying triplet insta- 
bilities, such as Bc2- Following the technique proposed 
by Furche [6], P,! r A can be obtained as 



pit- 



2\{A X -B X ) 1 ' 2 M- 1/2 (A X 



1/2 



,(15) 



with M A - (A x - B x y/\A X + B X )(A X - B A ) 1/2 and the 
singlet orbital rotation Hessians 

(A x )iajb — (e a — £i)8ijS a b 

+2\(aj\w 1 : e \tb)~\(aj\w 1 : e \te), (16) 

(B x )i a jb = 2X(ab\w^ e \ij)-X(ab\w^\ji), (17) 

where ti are the RSH orbital eigenvalues. 

This method will be referred to as RSH+RPAx. In 
the limit of n = 0, it reduces to the standard KS scheme, 
while for fi — > oo it becomes a full-range ACFDT RPAx 
approach (with HF orbitals). We note that at second- 
order in the interaction w^ e the RSH+RPAx reduces to 
the RSH+MP2 method [26, 31]. 

Computational details. Equations (14)-(17) have been 
implemented in the time-dependent DFT development 
module [32] of MOLPRO 2008 [33]. We perform a self- 
consistent RSH calculation with the short-range PBE xc 
functional of Ref. 24 and add the long-range RPAx cor- 
relation energy calculated with RSH orbitals. The range 
separation parameter is taken at fi — 0.5, in agreement 
with previous studies [28], without trying to fit it. The 
A-integration in Eq. (14) is done by a 7-point Gauss- 
Legendre quadrature [6]. We use large Dunning basis 
sets [34] and remove the basis set superposition error 
(BSSE) by the counterpoise method. The full-range RPA 
and RPAx calculations have been done with PBE [35] 
and HF orbitals, respectively. The computational cost of 
RSH+RPAx is essentially identical to that of full-range 
ACFDT RPA. 

Results and discussion. Figure 1 shows RPA, RPAx and 
RSH+RPAx total energies with respect to the basis size 
for Ne2- In contrast with full-range RPA and RPAx, a 
fast convergence is observed for RSH+RPAx, similar to 
that of standard KS calculations. This improvement is 
explained by the fact that short-range correlations are 
compactly described by the short-range density func- 
tional. The reduced basis dependence of RSH+RPAx 
also means smaller BSSE [38]. Another important point 
illustrated by Fig. 1 is that the large RPA overestimation 
of the total energy, apparent for a large enough basis, is 
remedied by the RSH+RPAx method thanks to a more 
accurate description of short-range correlations. 

Figure 2 shows the interaction energy curves of Be2 
and Ne 2 . The RPA (with PBE orbitals) fails badly: a 
large bump for Be2 and an almost completely repulsive 
curve for Ne 2 are observed. A spectacular improvement 
is obtained with the RSH+RPAx method, which gives 
physically correct curves, especially accurate at medium 
and large distances. It also improves over both MP2 and 
RPAx. It has been verified that the short-range LDA xc 
functional of Ref. 21 provides very similar RSH+RPAx 
interaction curves, indicating a low sensitivity of the 
method with respect to the short-range functional. 

The proposed RSH+RPAx method overcomes many 
limitations of the RPA or related approaches. The results 



4 




show that it has the potential to describe successfully Acknowledgements. We thank F. Furchc, J. F. Dobson, 
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